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ABSTRACT 

The  laboratory  type  cesium  beam  frequency  standard,  Cs  1 ,  of  RRL  (Radio  Re¬ 
search  Laboratory)  features  a  hexapole  magnet  focusing  system,  a  55  cm  Ramsey 
cavity  using  a  coaxial  line-to-waveguide  transducer  and  a  digital  servo 
system. 

The  length  of  the  Ramsey  cavity  is  relatively  short.  Therefore,  in  order  to 
secure  a  good  signal-to-noise  ratio,  beam  uniformity  and,  to  decrease  the 
pulling  effect  caused  by  the  wings  of  neighboring  field-dependent  transitions, 
a  calculation  of  the  beam  trajectory,  using  the  Runge-Kutta  method,  was 
carried  out.  In  the  calculation,  the  following  items  were  considered:  1)  a 
beam  collimator  having  a  finite  area  instead  of  a  point  source,  2)  dependence 
of  the  effective  dipole  moment  on  the  position  in  the  magnet,  3)  motion  of  the 
atoms  in  the  IjJ  direction. 

After  evaluating  the  equations,  the  technical  advantage  of  the  focusing  system 
using  the  hexapole  magnets  became  clear.  The  results  of  the  calculation  agreed 
fairly  well  with  that  of  the  experiment.  This  recent  optimization  of  the  beam 
optics  has  probably  improved  the  microwave  power  dependent  frequency  shift  in 
the  accuracy  evaluation  of  RRL  Cs  1 . 

INTRODUCTION 

The  laboratory  type  cesium  beam  frequency  standard,  Cs  1,  of  RRL  features  a 
hexapole  magnet  focusing  system,  a  55-cm  Ramsey  cavity  using  a  coaxial  line- 
to-waveguide  transducer  and  a  digital  servo  system1 12.3,4,5, 

A  number  of  accuracy  evaluation  experiments  have  been  made  and  the  total  un¬ 
certainty  of  those  evaluations  has  been  gradually  lowered  because  of  our  ef¬ 
forts  towrd  various  kinds  of  improvement.  After  realizing  a  slow  velocity  beam 
optics  using  hexapole  magnets  of  3  cm  in  length  with  a  collimator  of  0.5  mm  in 
diameter,  the  microwave  dependent  shift  was  much  improved  and  the  total  uncer¬ 
tainty  of  1.1  x  10“ 13  was  obtained6.  Then  RRL  started  reporting  the  data  to 
the  BIH.  There  were,  however,  many  problems  remaining  unsolved  and  the  pos¬ 
sible  solutions  to  them  could  not  be  obtained  clearly.  Therefore,  a  calcula¬ 
tion  of  beam  trajectory  for  the  optimization  of  the  beam  optics  was  carried 
out  and  the  beam  optics  and  the  electronic  circuits  of  the  RRL  Cs  1  were 
improved.  The  magnetic  field  effect,  the  pulling  by  the  neaighboring  transi¬ 
tions  and  the  power  dependent  shifts  were  studied  experimentally  and 
theoretically7.  The  total  uncertainty  of  0.7  x  10“ 13  was  achieved  as  expected 
in  April ,  1986 . 
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In  the  following,  a  calculation  of  beam  trajectory  for  the 
optimization  of  beam  optics  is  explained  and  the  results  of  calculation 
are  compared  with  those  of  experiment. 

CALCULATION 

The  magnetic  potential  in  the  region  between  the  pole-pieces  of 
hexapole  magnet  is  given  by8 

V=Crscos34»  ( 1  ) 


where  the  cylindrical  coordinate (r,  z)  is  adopted  and  z  coordinate  is 
madeto coincide  with  the  center  axis  of  the  hexapole  magnet.  The  magnetic 
field  in  the  r  and  the  ^  directions  can  be  written  as  gradient  of  Eq.  (1) 


Hr  =  -HmaX(r/r°  )2C0S3'1' 

(2) 

M),„(r/r„)!si~3*  (3, 

IHI=4PUH7  - *■  Hmax  .tr„„  <4> 

are  obtained,  where  r0  is  a  radius  of  the  inscribed  circle  of  hexapole 
magnet.  The  force  on  the  atom  has  the  components, 

Fr =“»[  d  2  r /dt*-r  (  dTl>/dt )  2  ] 

(5) 

= m  1  /  r  (  r  2  d  *  /  d  t ) 

(6) 

where  m  is  the  mass  of  atom.  The  interaction  energy  is  given  by  Eq.  (4) 

U^*lH=»ieff*liHl=Ueff*Ho.ax(r/ro  >2 

(7) 

where  JU  is  the  magnetic  dipole  moment.  The 
scalar  potential,  so  that 

interaction  energy  is  a 

F=-V(J=U  _  (r/re)2 

eff  3r"maxv 

(8) 

where  £  is  the  unit  vector  of  the  r  direction, 
be  written  as 

Eq .  ( 5 )  and  Eq .  ( 6 )  can 

m[d2r/dt2-r(dil)/dt)2 

(9) 

^(r2d4./dt)=0 

( 10) 

The  effective  dipole  moment  is  a  function  of  H  and  has  a  different  value 
for  each  state  of  cesium  atom: 

ueff(F=7/2±1/2.mF) 

"  2  (1  +  1/ 2m  i+x2)1-""2  ^  “*fa) 

r 

(11) 
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where  x=-2g0  *H/  W=-3.052  x  10- 4 H  and  H  is  the  magnetic  field  intensity.  It  is 
difficult  to  solve  eqs.  (9),  (10)  and  (11)  analytically.  Therefore  they  are 
solved  as  a  group  of  differential  equations  (r,  r,  IjJ,  lfl)  numerically  by  the 
Runge-Kutta  method.  For  example,  the  numerical  values  of  (r,  r,  IjJ,  ijj)  are  cal¬ 
culated  30  times  for  the  approximation  in  the  case  of  the  6  cm  magnet.  It  is 
assumed  that  the  flight  time  in  each  of  30  segments  does  not  vary  with  the  r 
position,  and  the  motion  of  the  atom  under  the  condition  r<l/10ro  is  assumed 
to  be  the  -same  as  that  in  the  free  potential  space,  since  the  magnetic  poten¬ 
tial  in  the  region  is  very  small.  Hence,  the  equation  can  avoid  divergence  un¬ 
der  the  condition  r<0. 

In  the  calculation  the  beam  collimator  has  finite  area  instead  of  a 
point  source  and  it  is  assumed  that  many  collimators  of  a  0.1  mm  diameter 
gather  to  form  a  bundle  and  the  directivity9  of  each  0.1  collimator  is  the 
same  as  that  of  the  original  one. 

OPTIMUM  DESIGN  OF  BEAM  OPTICS 

Figure  1  shows  a  velocity  dependence  of  the  cesium  beam  intensity  at  the 
exit  of  the  collimator.  The  intensity  ratio  of  the  100  m/s  beam  to  the  most 
probable  velocity  270  m/s  beam  is  1:6.2.  The  size  of  the  collimator  may  be 
chosen  so  that  the  compensation  for  beam  intensity  by  increase  of  temperature 
does  not  induce  poor  beam  directivity  of  the  collimator.  Figure  2  show  a  r2 
-dependence  of  the  magnetic  field  intensity  between  the  pole  tips  of  the  6  cm 
hexapole  magnet  which  was  use  during  the  initial  phase  of  the  experiment.  The 
inscribed  radius  is  1.5  mm.  In  the  calculation  the  effect  of  deviation  from 
the  r2  law  is  considered.  Figure  3  shows  the  typical  calculation  results  of 
velocity  distribution  in  beams  using  a  6  cm  magnet  and  a  3  cm  magnet.  One  pair 
of  3  cm  magnets  was  selected  for  an  optimum  focusing  system,  since  the 
microwave  power  dependent  frequency  shift  was  much  improved  with  3  cm  magnets 
in  May  1984.  The  collimators  of  0.5  mm  diameter  and  1  mm  length  were  used. 
However,  the  efficiency  of  the  beam  using  these  collimators  was  low 
experimentally.  In  the  design  of  the  beam  optics,  the  following  items  were 
considered:  1)  beam  efficiency  for  better  frequency  stability,  2)  beam 
uniformity,  which  is  also  influenced  by  the  gravitational  force  on  the  atoms, 
and  its  effect  has  a  relation  to  the  beam  asymmetry  in  the  interaction  region 
of  the  Ramsey  cavity  and  in  the  cutoff  waveguide,  3)  pulling  effect  by  the 
wings  of  the  neighboring  transition  (F=4,  mf=±l)  to  (F-3,  mf=±l). 

Figure  4  shows  the  distribution  of  arrival  of  the  atoms  in  the  r- 
direction  at  the  detector  for  different  collimator  diameters.  The  distance  be¬ 
tween  the  collimator  and  the  magnet  is  2.3  cm  and  the  diameter  of  the  beam 
stopper  is  1.2  mm.  In  the  case  of  flop-out  type  detector,  Figure  4  indicates 
that  the  radii  of  0.4  mm,  0.6  mm  and  0.75  mm  are  optimum  for  the  collimators 
of  0.5  mm,  0.8  mm  and  1.0  mm  in  diameter  respectively.  Figure  5  shows  the 
velocity  distributions  for  different  collimator  diameters.  The  dependence  upon 
the  distance  between  the  oven  and  the  collimator  is  also  calculated  and  the 
same  characteristic,  but  the  dependence  of  the 
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probable  velocities  upon  the  distance  is  stronger  than  that  upon  the 
collimator  diameter.  Figure  6  shows  the  examples  of  Ramsey  resonance  at 
3  dB  below  the  optimum  power  for  different  collimator  diameters.  In  the 
case  of  0.8  mm  collimator  diameter,  the  linewidth  and  the  peak-to-valley 
of  Ramsey  resonance  at  the  temperature  of  100  eC  are  105  Hz  and  2.26  pA, 
respectively.  The  total  number  of  atoms  which  arrive  at  detector  without 
transition  is  about  7.3  times  as  many  as  the  number  of  the  (F=4,  m^O) 
energy  state  for  the  three  different  collimators.  Figure  7  shows  the 
estimated  short  term  stability  and  most  probable  velocity.  On  the  basis 
of  above-mentioned  calculations  and  previous  experimental  results,  the 
beam  optic  parameters  and  dimensions  are  decided  as  shown  in  Figure  8. 

EXPERIMENTAL  RESULTS  AND  DISCUSSION 

Figure  9  shows  the  obsereved  Rabi  transitions  at  3  dB  below  the 
optimum  power.  The  magnetic  field  is  usually  set  for  1.088188  T  5  in 
operation.  Figure  10  shows  the  beam  current  characteristic  around  the 
clock  transition.  Figure  11  shows  the  Ramsey  resonance  of  clock 
transition.  Its  linewidth  and  peak-to-valley  beam  current  at  an  oven 
temperarure  of  about  112  0  C  are  100  Hz  and  4  pA,  respectively.  The  vapor 
pressure  of  10CTC,  112  'C  and  140 0 C  in  oven  are  0.55X10“ 3  mmHg,  0. 90X10' 
mmHg  and  2.5X10“ 3  mmHg,  respectively  .  Judging  from  these  data  and  the 
results  in  Figure  7,  the  beam  efficiency  has  been  improved  3~4  times, 
compared  with  the  previous  one. 

Figure  12  shows  the  velocity  distributions  of  three  transitions, 
which  indicates  the  reproducibilities  in  the  same  and  opposite  beam 
directions.  The  improvement  of  beam  stoppers  in  their  pattern  and 
position  adjustment  were  also  effective  in  these  reproducibilities. 
These  results  are  satisfactory  and  agree  fairly  well  with  the  calculated 
ones,  as  shown  in  Figure  5-  Figure  13  shows  the  recent  microwave  power 
dependent  frequency  shift  together  with  the  results  in  the  previous 
experiment.  The  uncertainty  of  0.7X10“  is  obtained  in  the  accuracy 
evaluation  experiment.  It  was  1.1X10”  in  the  previous  one. 

CONCLUSION 

The  authors  have  designed  the  optimum  beam  optics  on  the  basis  of  a 
caluculation  using  Runge-Kutta  method  and  the  previous  experimental  data. 
The  calculated  results  agree  fairly  well  with  those  of  the  experiment. 
The  obtained  improvements  in  the  short  term  stability  and  the  uncertainty 
of  RRL  Cs  1  are  satisfactory. 

As  to  the  Runge-Kutta  method,  the  authors  adopted  this  to  design  the  beam 
optics  using  6  cm  magnets  at  the  initial  stage  of  experiment  1  and  to 
design  the  present  beam  optics  using  3  cm  magnets  .  Huang  Bingying  of 
National  Institute  of  Metrology  in  China  also  developed  this  method  to 
design  a  new  cesium  standard  using  hexapole  magnets  and  concluded  that 
this  method  was  useful  for  the  simulation  of  beam  optics12.  In  our 
calculaton,  it  is  a  demerit  to  take  much  time  to  perform  a  numerical 
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calculation,  for  example  15^20  minutes  typically.  Therefore  further 

improvement  for  efficient  and  precise  calculation  will  be  desirable. 
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Fig.  1  Velocity  distribution  at  the  exit  of  collimator. 

Fig.  2  r  dependence  of  the  magnetic  field  intensity. 

Fig.  3  Typical  velocity  distributions  in  beams  using  3  cm  magnets  and  6 
cm  magnets 

Fig.  4  Distributions  of  the  arrival  atoms  in  the  radius  direction  at 
detector. 

Fig.  5  Velocity  distribution  for  different  collimator  diameters. 

Fig.  6  Ramsey  resonance  for  different  collimators. 

Fig.  7  Estimated  short  term  stability  and  most  probable  velocity. 

Fig.  8  Improved  beam  optics  and  the  structure  of  Ramsey  cavity. 

Fig.  9  Observed  Rabi  transitions. 

Fig.  10  Beam  current  arround  the  clock  transition. 

Fig.  11  Ramsey  resonance  of  the  clock  transition. 

Fig.  12  Velocity  distribution  of  three  transitions. 

Fig,  13  Microwave  power  dependent  characteristic  on  beam  reversal, 
compared  with  that  in  the  previous  optics. 
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Fig.  3  Typical  velocity  distributions 
in  beams  using  3  cm  magnets  and  6 
cm  magnets 
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Fig.  4  Distributions  of  the  arrival 
atoms  in  the  radius  direction  at 
detector. 
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Fig.  6  Ramsey  resonance  for  different 
collimators . 
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Fig.  7  Estimated  short  term  stability  and  most  probable  velocity. 
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Pi_g,  8  Improved  beam  optics  and  the  structure  of  Ramsey  cavity. 
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RELATIVE  EREQUENCY  [  KHz  ] 

Fig.  9  Observed  Rabi  transitions 
transition. 


Fig.  10  Beam  current  arround  the  clock 
transition. 


Fig.  11  Ramsey  resonance  of  the  clock 
transition. 


Fig.  12  Velocity  distribution  of  three 
transitions. 
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RELATIVE  MICROWAVE  POWER 
1.0  :  OPTIMUM  POWER 

Fig.  13  Microwave  power  dependent  characteristic  on  beam  reversal 
compared  with  that  in  the  previous  optics. 
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